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Abstract 

The A baryons with a single heavy flavor which transfer the quark po- 
larization, have been studied both theoretically and experimentally. The 
S's with two heavy constituents are well treated in quark-diquark model. 
In this work we study the production of triply heavy baryons in the per- 
turbative QCD regime and calculate the fragmentation functions for ^ccc 
and in the c and h quark fragmentation, respectively. We then ob- 
tain the total fragmentation probability and the average fragmentation 
parameter for each case. 
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1 Introduction 



The quark model of hadrons has proved to be successful in describ- 
ing hadrons and their properties. In the heavy quark sector it predicts 
hadrons having c, h and t quarks as constituents. However, the discovery 
of the top quark [1] and the determination of its lifetime [2] made it clear 
that it cannot participate in strong interactions and therefore only the c 
and h flavors are left to take part in the hadron production interplay. 



Recently meson states constituting heavy flavor have received consider- 
able attention. Specially Be and B* states with be quark content have 
been in focus both theoretically [3] and experimentally [4] in the last few 
years. It is established that the fragmentation functions describing their 
production mechanism are calculable in perturbative QCD [3] and hence 
the total fragmentation probabilities and production cross sections are 
calculated in due course. 

Baryon states with heavy flavor fall into three categ ones. States con- 
taining one heavy flavor such as Ac and Aj, are interesting states due to 
the fact that they carry the original heavy flavor polarization. They are 
presently being studied experimentally [5]. The second category involves 
baryons with two heavy flavor like the states Sec, '^hh and H^c [6]. They 
are treated within the approximate quark-diquark model [7]. The model 
treats the production of the so called diquark perturbatively similar to 
the states such as Be- Then, it can be proved that the formation of a 
baryon out of the diquark is almost the same as the fragmentation of an 
antiquark into a meson. In this way one obtains the fragmentation func- 
tions, the total production probabilities and other relevant parameters 
which specify their properties. In the third category, we have baryons 
with three heavy constituents. If we follow the scheme used in the case of 
heavy mesons and assume that their fragmentation functions are calcula- 
ble in the perturbative regime, then we can calculate Feynman diagrams 
like the one in flgure 1 to obtain the fragmentation functions. There are 
eight such diagrams in the lowest order contributing triply heavy baryons, 
i.e. Occc, ^ccb, ^cbb and Qbbb production [8]. 

In this paper our aim is to calculate the fragmentation of the ^ccc and 
baryons in the lowest order perturbative regime and obtain their 
fragmentation functions in an exact analytical form. 

2 Kinematics 

We consider the fragmentation of a heavy quark Q into a QQQ system 
with three identical flavor. This procedure is illustrated in Figure 1. We 
have used an inflnite momentum frame in which all of the particles are 



Fig. 1. Feynman diagram illustrating the lowest order fragmentation of a heavy quark, Q, into 
a ^QQQ baryon. The four momenta are labelled. 

moved in the forward direction, i.e. the longitudinal direction along the 
z axes, where the QQQ moves. We let the original quark keep its trans- 
verse momentum. Furthermore, we assume that the two antiquark jets 
move almost in the same direction. This assumption is justified due to 
the fact that the very high momentum of the initial heavy quark will 
predominantly be carried in the forward direction. Due to momentum 
conservation, the total transverse momentum of the two jets will be iden- 
tical to the transverse momentum of the initial quark. In this context the 
four momenta of the particles will assume the following form 



We have used the fragmentation parameter, z, as defined in the literature, 
i.e., 



(1) 



z = 



{E + p{)b _ Eb 



(2) 



(E+1>||)q Eq- 



The last step follows form application of the infinite momentum frame. 



Therefore, the final state particle energies are parameterized as follows 



Po = ZPo, 



(3) 



where Po = To + Qo + So is the energy of the baryon. Therefore 



XlZPo, 



Qo = 



X2ZP0, 



xszPo- 



(4) 



Here the x's are the energy ratios carried by the constituents. Since the 
constituents are identical and fly together, it is found that xi = X2 = 
= 1/3. This is consistent with our argument about the wave function 
for such states in the next section. We also have assumed that the two 
anti-quarks which initiate the two jets have equal energies, i.e. 



On the other hand due to our discussion about transverse momentum, 
we have 



We will discuss this later assumption in the final section. 

3 Calculation of the Fragmentation Functions for 0,qqq 

We are now ready to calculate the diagram shown in figure 1. The frag- 
mentation of a heavy quark Q into a heavy baryon ilggg is obtained by 
squaring the total amplitude and integrating over final state phase space. 



ko = K = -il- z)p'^. 



(5) 




(6) 



(7) 



where Tg is the amphtude of the baryon production which involves the 
hard scattering amphtude Th and the non-perturbative smearing of the 
bound state. The average over initial spin states and the sum over final 
spin states are performed. The heavy hadron production amplitude is 
composed of a partonic part, which can be calculated using perturbative 
QCD, and a non-perturbative part, which describes the transition of free 
quarks into the final state hadron. In the framework of non-relativistic 
quark model, this non-perturbative part could be accounted for through 
the wave function which is calculable using potential models. Since at 
present there is no known information concerning such wave functions, 
we have assumed a delta function type wave function for them. This 
assumption guarantees that the constituents will fly parallel and have 
no transverse momentum with respect to their direction of motion. This 
is also consistent with our assumptions in section 2. The hard scatter- 
ing amplitude which is obtained by perturbative calculations of the tree 
diagram in figure 1 , may be put in the following form [9] , 



Here as = g^/^'n is the strong interaction coupling constant and F in- 
dicates that part of the amplitude which embeds spinors and gamma 
matrices. The 1/p's are the propagators of the two gluons and the inter- 
mediate fermion respectively. 

To absorb the soft behavior of the bound state into hard scattering am- 
plitude we have used the scheme introduced in [7]. The probability ampli- 
tude at large momentum transfer factories into a convolution of the hard- 
scattering amplitude Th, and baryon-distribution amplitude 0m [10], i.e.. 



Th 



(8) 



TB{ki,Pi) = J[dx]THiki,pi,Xi)(t)B{xi,q''^) 



(9) 



where Th is given by (8) and 0s is the probability amplitude to find 
quarks co- linear up to a scale q'^ in the baryonic bound state. In (9), 
X-I s are the momentum fractions carried by the constituent quarks and 



[dx] = dxidx2dxsS{l — xi — X2 — x^). In view of our early discussion 
in this section, we propose the following expression for the probability 
amplitude 



= fB S{x, - (10) 

where is the baryon mass and /b refers to the characteristics of the 
baryon bound state and is similar to the meson bound state where the 
decay constant /m is introduced. Putting this expression and (8) in (9) 
and carrying out the necessary integrations, we find 

^ 247rV,m'fBCF T 

^ l2p'^p^k'^ko gi{z)g2{z)g3{z){p^ + ko-\-K-p'o)' 



Now we are able to obtain the fragmentation function in (7) as 
(487r2a>VBCi.)2 



Diz,fi) = 



8 

f \ Es rr^3(p + k + k' - pOdWkd^k' 

X / - — ^ 2- (12) 

PoV'okoK[gi{z)g2{z)gz{z){p^ + ko + k'^- p'J] 



Spin sum-average of TV for Figure 1 is easily calculated using the RE- 
DUCE. To do the phase space integrations in (12), first we consider the 
integral. 



/ = 



S^{p + k + k' -p')d^p 

P'oiPo + ko + k'^ - p'^y 

^° (13) 



where 

^ ^ 3m2 z 3 V 4m2/ 1- z ^ ^ 



Also we note that 



/ /(^,kT)A = / f{z,kT)dkLd'kT 



= m'kJiz, {kr')) = mXfiz, (15) 

and 

/ /(^,k'T)A' = / f{z,k'T)dk'Ld% 

= mXfiz, {k'^')) = mXfiz, ^{p'r')). (16) 

Here, instead of performing transverse momentum integrations, for sim- 
plicity we have replaced them by their average values. Putting all this 
back in (12), we obtain the fragmentation function as, 

4 4 f2 p2 



108m2^4(l - zYfizygizf 

X k^z^ + 4^^/(83 - 130z + 51^2) 

+6^^^(1413 - 3084z + 3022^2 - 2156z^ + 821z^) 
+4^2^2(18711 - 51678Z + 69417^^ - 70308^^ 
+53529^^ - 25950z^ + 6343z*^) + 222345 - 740664^ 
+1179036^2 _ 1253448^^ + 90126/ - 388872^ 



+109916^^ - 49912z^ + 20649z^ 



(17) 



Here we have defined ^ = {p%) /w?. g{z) comes from the propagators and 
have the following form 



3 4/ z 



and f{z) is due to the energy denominator given by (14). Replacement 
of f{z) by g{z) which we have done in the original manuscript, changes 
the fragmentation function only slightly. 



The fragmentation function Dc^^^^^{z, fio) and Di,^q^^^^^{z, fio) are easily 
obtained from the above by letting m = rric^mb and using appropriate 
/s, (y-s and jio values. 



4 Results and Discussion 

We were able to calculate the process of direct c and h quark fragmen- 
tation into Q^ccc and r2?,56 baryons. In doing so we had to follow certain 
assumptions. Firstly we have considered only the dominant contributing 
Feynman diagram in leading order. This assumption reduced the com- 
plexity and the length of the calculation and enabled us to obtain analytic 
forms of the fragmentation functions. Our second assumption concerns 
kinematics. We believe that the high momentum of the process has to be 
taken away in the forward direction and let the two antiquarks carry the 
transverse momentum of the initial heavy quark. Furthermore since they 
are identical, we have considered equal contribution from them both in 
magnitude an in direction. Therefore, we have established equations (5) 
and (6) and used them in our calculation. To see how our later assump- 
tion works, we have set the kinematics by allowing k and k' to share the 
jet energy-momentum. We have let k = x{l — z)p' and k' = {l — x){l — z)p' 
where x is a variable which is between zero and one. We have repeated 
our calculations and studied the behaviour of ilccc fragmentation function 
with the same parameters as before. It is revealed that as x increases, 
the function grows rapidly and gives the highest peak at x = 1/2. As x 
increases further, the peak falls rapidly. Since there is not much informa- 
tion about the wave functions of the triply heavy baryons at hand, we 
have reduced the non-perturbative smearing of the bound state to a delta 
function times a factor which is much like the meson decay constant. We 
have denoted this constant by /b and assumed to take 0.25 GeV both 
for Occc and baryons. 

In obtaining (17) we have not performed the transverse momentum inte- 
grations. Instead we have replaced the variables by their average values. 
However the numerical integration converges well for sufficiently large 
transverse momentum. Let us now sketch the behaviour of our frag- 
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Fig. 2. The behavior of J^ccc (sohd) and ri^bb (dashed) fragmentation function at the respective 
fragmentation scale. 



mentation functions. Figure 2 shows the behaviour of Dc^^^^^,{z^ fio) and 
Db-,nbbh{^^ A^o) in the fragmentation scale fio. In drawing them we have as- 
sumed that rric = 1.25 GeV and = 4.25 GeV. The scales are /io=6.25 
GeV and /Uo=21.25 GeV respectively. We have set )= 1 GeV and 
included the colour factor of Cf=7/6 obtained using color line count- 
ing rule. Consistent with the study of Be and B* states, we have taken 
as = 0.26 for Qccc and as = 0.18 for Qbbb [3]. 

At leading order in a, one has Jq Pq^q{z, fi)dz = [12], and the evolu- 
tion equation implies that the fragmentation probability Jq Dq^b{z, fi)dz 
does not evolve with the scale fi. Therefore, the fragmentation probabil- 
ity is a universal characteristic of the production rates. The evolution 
only moves the z-distribution to small values of z. We have obtained this 
quantity for Qccc and Qbbb using our fragmentation functions. The other 
relevant kinematical parameter is the average fragmentation parameter. 
Our results for the fragmentation probabilities and {z) appear in Ta- 
ble 1. It is seen that our analysis give very close {z) values for flccc and 



Tabic 1 

Fragmentation probability and {z) for different states. 



^ccc 


^bbb 


Prag. Prob. 2.789 x 10^ 


6.459 X 10'^ 


(■■) 0.522 


0.535 



Q^hbb- The fragmentation probabilities in Table 1 suggest that considerable 
event rate is expected both at the Tevatron and the LHC . 
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